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ABSTRACT 

The present study experimentally demonstrated that solid H2O is formed 
through the surface reaction OH + H2 at 10 K. This is the first experimental 
evidence of solid H2O formation using hydrogen in its molecular form at temper- 
atures as low as 10 K. We further found that H2O formation through the reaction 
OH + H2 is about one order of magnitude more effective than HDO formation 
through the reaction OH + D2. This significant isotope effect results from differ- 
ences in the effective mass of each reaction, indicating that the reactions proceed 
through quantum tunneling. 

Subject headings: astrobiology - astrochemistry - atomic processes - ISM: clouds 
- ISM: molecules 



1. Introduction 

1.1. H2O formation 

Water (H2O) is the most abundant solid component of icy grain mantles in molecular 
clouds. Previous theoretical studies demonstrated that the observed abunda nce of solid H2O 



in mo lecular clouds cannot be explained by gas-phase synthesis alone (e.g. iHasegawa et al. 



1992h. Therefore, H2O i s mainly considered to be formed on the surfaces of interstellar grains. 



Tielens fc HagenI (|1982[ ) proposed a chemical reaction network by which the formation of 
solid H2O begins with hydrogenation of O atoms, O2, or ozone, O3, on the grain surfaces. 
Since this proposal, many researchers have experimentally explored possible pathways for 
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the formation of solid H2O under the simulative conditions of dense molecular clouds. These 
studies showed that solid H2O is produced through various surface reactions at temperatures 
as low as 10 K. The main surface reactions leading to the formation of solid H2O and some 
related to H2O formation are listed in Tabled] 

The simplest pathway toward the formation of solid H2O is the sequential hydrogenation 
of O atoms: 



+ H OH, 
0H + H->H20. 



(1) 
(2) 



Cuppen fc HerbstI (120071 ) proposed that reaction ([2]) is the main route toward H2O formation 
in diffuse and translucent clouds in which H atoms are more prevalent than H2 molecules. 
Several groups have experimentally appro ached these reactions using D atoms instead of H 
atoms. For example, iHiraoka et al.l ( 119981 ) sprayed D atoms over O atoms trapped in an N2O 
matrix at 12 K. The O and D atoms were prod uced by dissociation of, respectively, N2O 
and D2 under a dc discharge. iDulieu et al.l (I2OIOI ) codeposited O and D atoms produced by 
dissociat ing O2 and D^, r espectively, onto non-porous amorphous H2O ice at 10 K. More 
recently, IJing et al.l (120 111 ) deposited D atoms and ^^0 (instead of ^^O) at 15-25 K on a 
synthesized olivine surface. All three groups claimed that D2O (HDO) formed during their 
experiments through surface reactions ([1]) and ([2]); however, other 0-related species, such as 
O2 and O3, were also present in or on the samples, and it was difficult to separate information 
relating to reactions ([T]) and ([2]) from the different competing H20-formation channels related 
to O2 and O3, shown below. Further studies are desired to isolate these reactions ([T]) and 
([2]) from reactions involving species other than H and O. 

The second simplest formation pathway involves the sequential hydrogenation of O2 via 
the formation of hydrogen peroxide (H2O2): 



02 + H^H02, 
HO2 + H ^ H2O2, 
H2O2 + H ^ H2O + OH. 



(3) 
(4) 
(5) 



Miyauchi et al.l (120081 ) first revealed experimentally that these reactions proceed at 10 K 
and that a significant isotope effect is apparent for reaction ([5]) but not for reactions ([3]) and 
This is consistent with the fact that reactions ([3]) and (jlj) have essentially no barrier 



(IWalch et al.l Il988l: ISellevag et al.ll2008l ). whereas reaction ^ has a significant barrier of 
about 2000 K (IKoussa et al.ll2006l ) and requires quantum tunneling (IMiyauchi et al.ll2008l ). 



After the experimental study by iMiyauchi et al.l (120081). extensive stu dies were per- 
formed on the hydrogenation of O2 at low temperatures. Iloppolo et al.l (120081 ) revealed that 
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reactions (I3])-(I5]) occur even at 28 K. Furthermore, it was found that the structure of the 



H2O ice formed through reactions ([3]) -(151) is amorphous (lOba et aLll2009l ). which is consistent 
with the astronomical observations (ISmith et al.lll989[ ). The recent detection of H2O2 in the 



cloud core of p Op h A further supporte d the occurrence of these reactions on the surface of 
interstellar grains ( Bergman et al. 12011) because gas-phase synthesis of H2O2 through reac- 
tions dni) and dlD is not favored (IMousavipour fc Sahebll2007l ). In the gas phase, H2O2 formed 



through reaction (1^ further dissociates into two OH species using the excess energy from 
H2O2 formation (IMousavipour fc Sahebl 120071 ). In the solid phase, the release of the excess 
energy from reaction ^ into the cold surface will stabilize the formed H2O2 on the surface. 
Even though the dissociation occurs, OH radicals are expected to recombine, resulting in 
the recovery of H2O2 (reaction dH})) and/or the formation of H2O (reaction d?])), as follows: 



OH + OH^H202, 
0H + 0H^H20 + 0. 



(6) 
(7) 



The reaction HO2 + H in the gas phase was theoretically predicted to follow two channels 
to produce H2O via an intermediate species as follows: 



HO2 + H ^ H2OO* ^ H2O + 0(^D), 
HO2 + H ^ H2O..O ^ H2O + 0(^P), 



(9) 



wher e H2OO* and H2O..O denote short-lived intermediate species ( IMousavipour &: Saheb 
20071 ). Reaction channel dHD as a whole is slightly exothermic by 10 kJ mol~^, whereas 
the latter part of the reaction (H^OO * — > H2O + 0(^D)) is highly endothermic by 151 
kJ mol~^ (IMousavipour &: Sahebl 120071 ). indicating that reaction channel ^ would not go 
to completion if the intermediate H2OO* dissipated the excess energy to a cold surface. 
Reaction has a huge barrier of about 9300 K, implying that the reaction proceeds only 
with difficulty, even with tunneling under astrophysically relevant conditions. Furthermore, 
the branching ratio to pr oduce O atom i s only 2% for reaction of HO2 + H in the gas phase 
(lKevserlll986l ). Although ICuppen et al.l ( 120 lOl ) claimed that the reaction HO2 + H ^ H2O 
-|- O occurred in their O2/H codeposition experiments, this reaction would be unlikely in the 
present experiment. 



Reactions of HO2 with H may also lead to the formation of H2 and O2 ('^S^ ) ( IMousavipour fc Saheb 



20071 ): 



H02 + H^H2 + 02('Sr). 



(10) 



In the gas phase, this reaction occurs much less frequently than reaction dlj) (lKeyserlll986l ). 
Even if O2 were reproduced through reaction (ITOl) . it must be hydrogenated once again to 
yield HO2 and finally H2O2 like in the 02-hydrogenation experiments. 



-4- 



The sequential hydrogenation of ozone (O3) may occur in molecular clouds: 



O3 + H -> O2 + OH. 



The products (O2 and OH) may be further used for H2O formation. iMokrane et al.l (120091 ) 
investigated these reactions at 10 K by irradiating solid O3 on H2O ice with D atoms. The 
react ion products wer e anal yzed using temperature-programmed desorption (TPD) meth- 
ods. iRomanzin et al.l (I2OIII ) performed experiments of the hydrogenation of O3 at higher 
temperatures (25, 40, and 50 K), and the reaction products were monitored in-situ by Fourier- 
transform infrared (FTIR) spectroscopy. Both studies provided results indicating that the 
formation of H2O was initiated by reaction (fTTi) . 

Note that OH is always produced along each of these three pathways: reaction ([T]) for 
the hydrogenation of O atoms, reaction for that of O2, and reaction (ITT!) for that of O3. 
Hence, in addition to the H-atom related reactions shown above, the following H2-related 
reaction would lead to H2O formation: 



OH + H2 ^ H2O + H, 



(12) 



which has an activation barrier of ~2100 K in the gas phase ( lAtkinson et al.ll2004j ). Despite 
the large activation barrier of reaction ( 1T2|) . reaction f|T2|) is often thought of as the most 
common route to the formation of H2O on grain surfaces in dense molecular clouds in which 
typical temperatures are abo ut 10 K, the UV flux is very low, and H^ is the dominant H - 
bearing species over H atoms (iTielens &: Hagerull982l : lTielensll2005l : ICuppen &: Herbstll2007l ). 
At temperatures as low as 10 K, Arrhenius-type reactions having an activation barrier of 
~2000 K occur only rarely. Therefore, reaction (fT2l) requires quantum tunneli ng to proceed 



on gr ain surfaces, as observed in the gas-phase reaction at low temperatures (INguyen et al. 



2OIII ). The tunneling reaction is described in further detail in the next section. 



Reaction ([T2l) wa s experimentally studied in the gas phase (IRavishankara et al 
Oldenborg et al.lll992[ lAlagia et al.lll993l : iTalukdar et "al]ll996 ; iKrasnoperov fc Michaell l 



1981 



2004 



Orkin et al. 

mm- 

In most previous studies, OH was produced by photolysis of H2O or H2O2, 
or indirectly by photolysis of N2O to generate 0(^D) atoms, which further reacted with H2O 



(IVaghjiani et al. 



19891 : lOldenborg et al.lll992l : ITalukdar et al.lll996t lOrkin et al.l 120061 ). OH 



generated through those energetic processes should have signiflcant excess energy that should 
be sufficient to overcome the activation barrier of reaction (fT2|) . Hence, to investigate whether 
reaction (IT^ occurs through quantum tunneling on cold grains, non-energetic OH must be 
used for experiments. In recent experimental studies of low-temperature surface reactions, 
OH radicals were produced through various reactions and were further used in secondary 
reactions. For example, hydrogenation of O3 yielded OH as a reaction product (reaction 
ffTTj) ): however, OH formed in reaction ffTTj) could be highly energetic because the heat of the 
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reaction (~39,000 K; lKeyserlll979l ) is partitioned into OH (IRomanzin et al.ll201ll ). which is 
not suitable for investigating quantum tunnehng reactions. 

We recently succeeded in producing cold OH by dissociating H2O in a microwa ve-induced 



plasm a, follow ed by cooling to 100 K before reaching the reaction substrate ( lOba et al. 



2010al lbl. boilh . A similar technique for producing OH in the electronic and vibrational 
groun d states was applied to crossed beam studies of reaction ( fT2l) in the gas phase (lAlagia et al 
IQQsl ). where OH was produced in a radiofrequency-induced H2O plasma. We fou nd that the 



cold O H can react with CO and HOC O on a substrate even at 10 K to yield CO2 (lOba et al. 



2010al ) and H2CO3 (lOba et al.ll2010bl ). respectively. Moreover, it was experimentally demon- 
strated that cold OH reacts with another OH to y ield H2O2 and H2O through reactions 
and d?]), respectively, even at 40 K ( lOba et al.ll201ll ). Consequently, our experimental method 
is useful for studying chemical reactions related to cold OH on the surfaces of interstellar 
grain mantles. 

In the present study, we first performed the codeposition of cold OH with H2 on a cold 
substrate to study reaction (IT^ as a possible route to the formation of H2O on the surfaces 
of interstellar grains in dense molecular clouds. Next, OH and/or H2 isotopologues, such 
as OD and D2, were studied experimentally in an investigation of the isotope effect on the 
reactions to decisively conclude that these reactions occur via quantum tunneling. 



1.2. Quantum tunneling reaction 



The thermally activated reaction fll2p is not expected to occur on grain surfaces at 10 
K for the reasons given above. Previous experimental and theoretical studies proposed that 
quantum tunneling plays a significant role in reaction (fT2|) in the ga s phase, particularly 



at low temperatures (e.g. iTalukdar et al.l Il996t iNguyen et al.ll201ll ). IWatanabe fc Kouchi 



( 120081 ) described the importance of quantum tunneling in various chemical reactions on low- 
temperature grain surfaces in molecular clouds. Under the simple assumption that a reaction 
occurs through a rectangular barrier of height E^^ and width a on a potential energy surface, 
the rate of a quantum tunneling reaction (kq) is approximately represented by the following 
equation: 



(13) 



where vq and m denote the frequency of harmonic motion and particle mass, respectively. 
In the case of a bimolecular atom-transfer reaction represented by A + XB — t- AX + B, the 
center of mass translated from the initial state (A + XB) to the final state (AX + B) of 
a reaction on the potential energy surface is represented as the effective mass rric which is 
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defined by the following equation: 

mamb(l + c)^ + m\^m^c^ + mjn^ 



rrir 



M(l + c2) 



(14) 



where ma, m^, and mb are the masses of A, X, and B atoms, respectively, M is the mass 
of the linear triatomic molecule AXB (M = ma + rrix + mb), and c is the ra tio of the 
infini tesimal change of one bond distance to that of the other bond distance ([Johnston 
19661 ). When bimolecular atom-transfer reactions A + XB AX + B proceed through 
quantum tunn eling, an effective rn ass rric should be used to indicate the tunneling mass in 
Equation ([13]) jHidaka et ahlbnOQl ). 



Observation of an isotope effect on a reaction can be an indicator of quantum tun- 
neling. This is because tunneling effects strongly depend on the particle mass, such that 
a tunneling reaction involving a species with a lighter effective mass proceeds faster than 
with a heavier effective mass. Previous studies have investi gated the isot o pe eff ect on other 
astrophysically relevant chemical reactions around 10 K. iHidaka et al.l ([20071 ) performed 
hydrogenation/deuteration of solid CO at 10-20 K and found that hydrogenation of CO 
proceeds faster than deuteration of CO by a factor of about 12. H-D isotope exchange 
reac tions of methanol (C H.sOH) were also demonstrated to occur through quantum tunnel- 
ing ([Nagaoka et al.l 120071 ) . which would contribute to the observed abundance of deuterated 
methanol, such as CHD2OH. 

The importance of the effective rnass o n quantum tunneling reactions has been demon- 
strated experimentally. iHidaka et al.l (|2009[ ) extensively investigated H-D exchange reactions 
of formaldehyde (H2CO) at 10-20 K. The effective mass of the D-atom abstraction by H 
atoms from D2CO is ~1, whereas that of the H-atom abstraction by D atoms from H2CO is 
~0.5, which affects the rate of each reaction (IHidaka et al.ll2009t ). 



Most quantum tunneling surface reactions previously studied were H- and D-atom- 
related reactions. Previous experimental studies of H2-related quantum tunneling reactions 
have investigated the isotope effect on the reactions of methyl radicals (CH3) o r its deuterated 
species (e.g. CHD? and C D.s) in solid parahydrogen {Ea ~ 5300 K) at 5 K (IMomose et al. 



19981 : iHoshina et al.ll2004l ). Although these pioneering studies are important for understand- 
ing quantum tunneling reactions related to H2, their relevance to astrophysics is rather weak 
because their experiments were conducted under unfeasible conditions (in a parahydrogen 
matrix at 5 K) in molecular clouds. 

In the present study, we experimentally examined surface reactions related to H2 and 
OH isotopologues, as shown below, to investigate the isotope effects on reaction ([T2[) : 



OD + H2 ^ HDO + H, 



(15) 
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OH + D2 ^ HDO + D, 
OD + D2^D20 + D, 

OH + HD ^ H2O + D, 
OH + HD ^ HDO + H, 
OD + HD ^ HDO + D, 

OD + HD ^ D2O + H. 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



The use of isotopically substituted hydrogen (HD or D2) or hydroxy! radical (OD) in place of 
H2 or OH, respectively, can alter the effective mass of a reaction relative to that of reaction 
(12). In this case, the rate of a quantum tunneling reaction should differ from that of reaction 
(12) as well. Reported values of the activation barrier an d the calculat e d effe ctive mass for 
each reaction are summarized in Table [2l According to iHidaka et al.l ( l2009l ) , the value of 
c in Equation (HM was set to -1 in the present reaction system. An observation that the 
efficiency of each reaction varied according to the effective mass would strongly indicate that 
these reactions proceed through quantum tunneling. 



2. EXPERIMENTAL 



2.1. Apparatus and experimental conditions for the surface reactions 



Experiments were performed using the Apparatus for SUrface Reaction in Astrophysics 
system. Details of this apparatus and experimental conditions have been described pre- 



Naeaoka et al. 


2007; 


Oba et al. 


2010a. 


2011) 



radicals, together with H atoms, were produced b y the dissociation of H2O in a microwave 



discharge plasma (ITimmermans et al.l Il998l . Il999l ) in a Pyrex tube. H2O dissociation can 



lead directly to the formation of H2 and O. The species produced in the OH source are 
denoted "H2O fragments" in this paper. H2O fragments were transferred via a series of 
poly(tetrafluoroethylene) and a cold aluminum (Al) tube, which was connected to a He re- 
frigerator, and they impinged on the substrate. Initially, OH radicals may be excited when 
formed in the plasma, but after many collisions with the inner wall of a cold Al pipe, they are 
thermalized to 100 K. That is, they should be vibrationally and electronically in the ground 



(Alagia et al. 


1993; 


Kohno et al. 


2011) 



the next section. The deposition rate of OH radicals was not measured directly; the upper 
limit of the deposition rate should be 1.5 x 10^'^ radicals cm~^ s~^, which is equal to the 
deposition rate of H2O when the microwave source is turned off. H2 molecules delivered from 
a separate gas line through a capillary plate were codeposited with the H2O fragments onto 
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the substrate over 60 minutes. The flux of H2 was calculated to be 2.0 x 10^"^ molecules 
cm~^ s~^ based on the measured pressure inside the main chamber when filled with only H2. 

The reaction products were monitored in-situ by reflection-absorption-type FTIR with 
a resolution of 4 cm~^ in the spectral range 4000-800 cm~^. The desorbed species from the 
substrate were monitored using a quadrupole mass spectrometer (QMS). The temperature 
of the substrate, which was connected to another He refrigerator, was kept at 10 K in each 
experiment. The reactions in Table [2] did not occur if the codeposition experiments were 
performed at temperatures above 20 K. 

After codeposition, the TPD spectra were obtained by the QMS at a heating rate of 4 
K minutes"^. 

In the experiments using H2- and OH-isotopologues, the HD and D2 fluxes were the 
same as the H2 flux. QMS measurements and a resonance-enhanced multiphoton ionization 
(REMPI) method (see Section 2.2) confirmed that almost all (>95%) D2O or H2O introduced 
into the radical source was fragmented when the microwave source was turned on. On this 
basis, the deposition rate of OD was considered to be similar to that of OH because the flow 
rates of D2O and H2O into the radical source were adjusted to the same value. 



2.2. In Situ Detection of OH Radicals 



The direct detection of OH from the source was performed using anot her apparatus, 
named RASCAL, in our laboratory. RASCAL is described in detail elsewhere (IWatanabe et al 



2OIOI : iHama et al.ll201ll ). OH from the radical source was selectively ionized by (2 -|- 1) 



REMPI approximately 5 mm from the exit of the aluminum pipe and detected using a 
time-of-flight mass spectrometer. 

Figure [D^a) shows the REMPI spectrum of OH(v=0) via the D^J:-{v'=0)^ X^U{v=0) 
transition under conditions in which the Al pipe was equilibrated at room temperature. The 
presence of vibrationally excited OH(w=l) was sought using the RE MPI, with 3^E~(y ^ =0)-f- 



X^H(f=l) transition; however, no such evidence was obtained. iGreenslade et al.l (120051 ) 



calculated that the two-photon cross-sections of the D^S~(f' = 0) X^H(v=0) and 
3^S^(v'= 0) ^ X^H(f=l) transitions are comparable; therefore, the present results in- 
dicate that OH from the radical source was electronically and vibrationally in the ground 
state. The REMPI spectrum of OH(f =0) was assigned and modeled using the PGOPHER 
simulation program to determine the rotational structure (C. M. Western, available from 
|http:/ /pgopher. chm.bris.ac.uk). The const ants associated with the relevant electronic states 



were compared with those in the literature (jPieke fc Crosswhitdll962uGreenslade et al.ll2005 
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Huber fc Herzberglll979l ). Figure [U^b) shows the results of a spectral simulation of OH(f=0) 
assuming a rotational temperature of 300 K. By comparison with the spectral simulation, 
the OH(f =0) from the radical source was expected to be almost thermally equilibrated with 
the Al pipe before reacting on the substrate. 

In an additional step, the adsorbates were photodesorbed using the weak second har- 
monic radiation (532 nm) of a Nd:YAG laser and were analyzed by the REMPI method. We 
detected OH(f =0) signals in the species photodesorbed from the cold substrate (10-40 K) 
during the deposition of H2O fragments. Photolysis of H2O or H2O2 on the substrate, which 



could be potent ial sources for OH(i;=0), was negligible at 532 nm ( IChu &: Anastasid 12005 



Kobayashilll983l ): therefore, this result clearly indicates the presence of OH in the electronic 
and vibrational ground states on the substrate at 10 K. The rotational temperature of the 
photodesorbed OH(f =0) was about 100 K. This value is slightly higher than the substrate 
temperature. The OH radical is considered to gain a small amount of energy during photo- 
stimulated desorption. The present results are consistent w ith previous experimental results 
obtained using a similar method at 3.5 K ( IZins et al.ll201ll ). 



3. RESULTS AND DISCUSSION 
3.1. OH + H2 

The deposition only of H2O fragments on a substrate at 10 K resulted in the formation 
of H2O, H2O2, and O3, as identified in the product IR spectrum at 1650, 1400, and 1040 
cm~^, respectively (Figure [2](a)). Hereafter, this experiment is denoted as the H20-fragment 
depositio n. In contrast, t hese molecules were not observed in the H20-fragment deposition 



at 60 K (jOba et al.l 120111 ). If H2O, H2O2, and O3 formed in the gas phase, they should be 
adsorbed on the substrate even at 60 K, which is well below their desorption temperatures. 
We therefore concluded that these molecules were formed through surface react ions at 10 



K. T he formation pathways of these molecules are described in our recent paper (jOba et al. 



2OIII ). The peak area for each reaction product increased linearly with the deposition time. 
In addition to the three types of molecules listed above, QMS measurements revealed that 
H2 and O2, which are infrared-inactive, are present in the gases from the source and/or are 
produced on the surface. 

Figure [2](b) shows the IR spectrum of the reaction product obtained after codeposition of 
H2O fragments with H2. This spectrum was obtained under the same experimental conditions 
as those for the H20-fragment deposition: the deposition rate of OH, the experimental time, 
and substrate temperatures. The resulting spectrum was similar to that obtained in the 
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H20-fragment deposition, although the peak area associated with the OH bending of H2O 
(1650 cm~^) increased by a factor of two compared to that observed in the H20-fragment 
deposition (Figure |2|^c)). This result indicates that the introduced H2 was used for the 
formation of H2O through surface reactions with 0-related species, such as OH, O, O2, 
and/or O3. We concluded that H2O was formed by reaction of H2 with OH freaction (fT2|)') 



on th e surface at 10 K, despite the relatively large barrier E^, of about 2100 K ( lAtkinson et al. 
2004 ). for the following reasons. 



First, we confirmed that codeposition of H2 with cold O atoms produced by the dissoci- 
ation of O2 in a microwave-induced plasma did not result in a reaction at 10 K, except that 
O3 formed (Y. Oba et al., unpublished data). This result indicates that cold O atoms in an 
H2O plasma will not react with H2 either. This is in good agreement with the fact that the 
reaction. 



O + H2 ^ OH + H, 



(22) 



i s endothermic by 8 kJ mol 



( iBaulch et a 



action (22) ( 



[= 960 K) with a significant activation barrier of 3160 K 
1992!) ■ although theore tical studies assumed H2O formation begins with re- 



Cazaux et al.ll2010l . 1201 ll ). Second, O2 never reacts with H2 under the present 



experiment al conditions. Third, O.s does not react with H2, as experimentally demonstrated 
previously (IMokrane et al.l l2009l : iRomanzin et al.ll201ll ). In addition, if HO2 is produced 
through reactions of the H2O fragments, it may react with H2 to yield H2O2: 



HO2 + H2 ^ H2O2 + H. 



(23) 



In this case, the H2O2 produced may further yield H2O through reaction ([5]). This, however, 
is unlikely in the present experiment because reaction f l23p is highly endothermic with a 
considerable activation barrier (Tabled]). 

The column density of the H2O product was calculated from the peak area of the 
OH bending band at 1650 cm~^ and the previously reported integrated band strengths, as 
described in iHidaka et al.l (120071 ). T he band strength used was 1.2 x 10^^'' cm molecules"^ 
which is by a transmission method (iGerakines et al.lll995l ). We confirmed that the ratio of 
the band strengths for OH stretching and bending bands of H2O is almost independent of 
the method (transmission or reflection). While, the difference in the absolute value of band 
strengths for the reflection and th e transmission meth ods was found to be within the factor 
of two in our experimental setup (IHidaka et al.l 120091 ). Moreover, it is difficult to determine 
the exact column density even by the transmission method in the layered and the mixed 
samples. However, since our discussion is based on relative reaction yields, the difference in 
the absolute value of the band strengths between the two methods little alters the present 
experimental results and discussion. 
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By subtracting the H2O band obtained in the H20-fragment deposition from that ob- 
tained in the codeposition experiment, we determined the increase in the column density of 
H2O through reaction (12) to be 3.5 x 10^^ molecules cm~^. Note that it is not obvious 
whether H2O yield obtained by the subtraction gives the total yield of H2O by the reaction 
OH + H2. In the H20-fragment deposition, two competing reactions produce H2O: reactions 
([7]) and (fT2|) . In the codeposition experiment, additional H2 molecules enhance reaction (fT2|) . 
That is, a part of OH amount which was used by reaction ((Tj) in the H20-fragment depo- 
sition may be consumed by reaction f ll2l) . Therefore, H2O column density obtained by the 
subtraction may be underestimated compared to the total yield of H2O by reaction f|T2|) in 
the codeposition. Further details on this issue are shown in the next section. 

Accompanied by the increase in the H2O column density, the peak area of the OH 
bending in H2O2 at 1400 cm~^ decreased relative to the peak area measured in the H2O- 
fragment deposition (Figure [2](c)). This is probably because the abundant H2 which was 
codeposited with H2O fragments inhibited the formation of H2O2 through reaction ([6]) (OH 
+ OH — ^ H2O2); the OH fragments instead formed H2O via reaction (fT2|) . Reaction ([2]) 
(OH + H — )■ H2O) was assumed to contribute to H2O formation only to a minor extent in 
the present study because most H atoms in the H2O fragments must be consumed by H-H 
recombination on the surface or few H atoms were present in the gases from the source. 
In fact, the codeposition of O2 with the H2O fragments did not increase the amount of 
H2O2 relative to the H20-fragment deposition, despite the barrierless mechanism for O2 
hydrogenation. 

The absorption at 1040 cm~^ due to O3 was not observed in the codeposition experiment 
(Figure [2]^b)). It is possible that O3 was consumed by H atoms (reaction (fTT!) ) that formed 
through reaction (fT2l) and/or the O atoms were similarly consumed by H atoms (reaction 
([T])) prior to the formation of O3 (O2 + 0—7- O3). Reactions ([T]) and (ITT!) yield energetic 
OH, which may further react with H2 to yield H2O through reaction (1121) . Even if this 
were the case, the low O3 column density produced in the H20-fragment deposition (0.6 x 
10^^ molecules cm~^) suggested that such reactions did not significantly contribute to the 
observed H2O column density. Nevertheless, a lack of O3 could provide indirect evidence that 
the reaction OH + H2 occurred. In fact, the codeposition of D atoms with H2O fragments 
did not yield O3 formation (see Section 3.3). 

3.2. OD + H2 

Similar experiments were performed using OD instead of OH. If OD reacts with H2, the 
formation of HDO is expected to occur through the following H-atom abstraction from H2 
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by OD: 



OD + Ho ^ HDO + H, 



(15) 



which has almost the same activation barrier and effective mass as reaction (lT2i) (INguyen et al, 
201 ll . Table 2). In the OH/H2 codeposition experiment described in the previous section, 
the reaction products were basically the same as those of the H20-fragment deposition ex- 
cept for the presence of O3. In contrast, HDO, the reaction product of reaction (IT^ . is not 
produced unless both H and D-related species are used, that is, HDO is never produced in 
the deposition only of H2O or D2O fragments. The formation of HDO should thus provide 
direct evidence that reaction f|T5|) occurs at 10 K. 

Formation of D2O, D2O2, and O3 was observed in the product IR spectrum upon the sole 
deposition at 10 K of D2O fragments (Figure [31(a)) (hereafter, this experiment is denoted 
as the D20-fragment deposition), which were produced in a manner similar to that used 
to produce H2O fragments. These results indicated that the isotopic substitution of H2O 
fragments did not significantly alter the reaction pathways on the cold substrate. 

The codeposition of D2O fragments with H2 resulted in the appearance of new peaks at 
3392, 2858, and 1453 cm~^ (Figure [3](b)). Absorption bands at 3392 and 1453 cm~^ were 
typically observed for solid HDO produced by depositing an H2O and D2O niixture on a 
substrate at 100 K or lower (IHornig et al.lll958l : iDevlin et al.lll986l : iGalvez et al.ll2011l). The 
absorption band at 2858 cm~^ may correspond to the 2^2 vibration of HDO (jDevlin et al. 



19861 ). Because the only H source is H2 and other 0-related species, such as O atoms and 
O3, hardly react with H2, as explained in Section 3.1, we conclude that H2 reacted with OD 
to yield HDO on the surface at 10 K. The presence of HDO in the reaction product was also 
confirmed by the TPD experimental results. Figure H] shows the integrated TPD yields of 
m/ z = 19 (HDO) plotted with the peak area of the band at 1453 cm~^. 

The column density of HDO was calculated from the peak area of the OH-stretching 
band, and the previously reported integrated band strength of 1.3 x 10~^^ cm molecule"^ 
for the OH-stretching band of the solid HDO ( llkawa fc Maedalll968l ) was found to be 4.7 x 
10^^ molecules cm~^. 



In contrast to the experimental result shown in the previous section, the obtained column 
density of HDO should represent the total yield of HDO by the reaction 0D+ H2 because 
HDO is not formed in the D20-fragment deposition. The HDO yield is very similar to the 
estimated yield of H2O by the reaction OH + H2 (3.5 x 10^^ molecules cm~^) in the previous 
section. It is reasonably considered that our estimated H2O yield represents the total yield 
of H2O by the reaction OH + H2. This is based on an assumption that reaction f lT2|) has the 
same efficiency with reaction ffTS]) (see Section 3.6 for further details). 
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3.3. OH + D2 

We next codeposited H2O fragments with D2 to investigate the efficiency o f reaction (fT6l). 



which had almost the same barrier as reactions f|T2l) and f|T5|) in the gas phase (iNguyen et al. 



201ll ). The effective mass of reaction f|T6|) was twice that of reactions (|T2|) and (|T5l) (Table 
Figure E] shows an IR spectrum of the reaction product obtained after codeposition of 
H2O fragments with D2. The spectrum resembled that obtained from the H20-fragment 
deposition (Figure [2]^a)); however, a broad peak grew in at 2475 cm~^, which is a typical 
position for OD-stretching bands. In ad dition, a small peak was observed at 2150 cm~^. 



This peak may be associated with D2O2 (ILannon et al.lll97ll ). which is probably produced 



through deuteration of O2 by the D atoms formed through reaction 0161) . Because the only 
D source is D2 in this experiment, the appearance of an OD-stretching band indicates that 
reaction f|T6|) occurred at 10 K. The presence of HDO was verified in the TPD experiment. 



The column density of the formed HDO was calculated based on the peak area of the 
OD stretching band at 2475 cm~^ and the previously reported integrated b and strengths of 



5.0 X 10 cm molecule for the OD stretching band of the solid HDO (llkawa fc Maeda 



19681 ). The column density was found to be 2.7 x 10^^ molecules cm~^, much lower than 
the column density of the H2O formed via reaction (fT2l) under the similar experimental 
conditions (see Section 3.1). A low efficiency was also inferred from the presence of O3 in 
the reaction product (Figure ED. As described above, if reaction fITBl) proceeds efficiently to 
produce significant quantities of D atoms, O3 would be consumed. In fact, the codeposition 
of H2O fragments with D atoms did not result in observation of O3. 



3.4. OD + D2 

Reactions of OD with D2 are expected to yield D2O and D atoms as follows: 



OD + Do ^DoO + D. 



(17) 



Because the effective mass for reaction (IT71) exceeds those of reactions (IT^ and (IT^ and 
is comparable to that of reaction f|T6|) . reaction f|T71) is expected to be less efficient than 
reactions f|T2|) and f|T5|) . Figure [6] shows an IR spectrum obtained after codeposition of D2O 
fragments with D2. The spectrum resembled that of the D20-fragment deposition (Figure 
IHH^a)). A slight increase in the OD bending band of D2O was observed at 1200 cm~^, which 
corresponds to an increase in the D2O column densi ty by 3.5 x 10^'^ mol ecules cm~^. The 
band strength used was 8.0 x 10~ cm molecule"^ (IMivauchi et all 120081 ). 
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3.5. OH or OD + HD 

Unlike reactions that involve H2 or D2, two reaction channels are available for OH in 
the presence of HD: 

OH + HD -> H2O + D, (18) 
OH + HD^HDO + H. (19) 



These reactions are associated with nearly identical barrier heights (iNguyen et al.l |2011[ ) 
whereas the effective masses differed considerably from one another (0.48 for reaction flTSl) 
and 0.90 for reaction (fT9|) ). The experimental results shown in Sections 3.1-3.4, therefore, 
predict that reaction (fT8|) is more efficient than reaction (fT9|) . The branching ratio for these 
reactions in the gas phase was experimentally estimated, suggesting that reaction (fT8|) was 



3-4 times more efficient than reaction fll9p (ITalukdar et al.l 119961 ). H2O fragments were 



codeposited with HD to investigate the efficiencies of reactions (1T5]1 and (IT^ . 

The formation of H2O and H2O2 was observed in the product IR spectrum obtained 
after codeposition of H2O fragments with HD. O3 was not observed (Figure [7|(a)). The peak 
area for H2O increased by a factor of 1.5 compared to the H20-fragment deposition (Figure 
[Tl^b)), which corresponded to an increase in the H2O column density by 3.1 x 10^^ molecules 
cm~^. In addition to the H2O increase, small peaks were observed at 2435 and 2114 cm~^ 
(Figure [7](b)). The former peak corresponded to an OD-stretching band, implying that HDO 
was formed through the reaction OH + HD; however, the observed OD-stretching band was 
expected to include large contributions from species other than HDO. If reaction fll8p is 
the dominant channel relative to reaction (IT^ . D atoms are produced together with H2O, 
which results in the formation of D2O2 through reactions between D and O2 contained in 
the H2 O fragments. In fac t, the peak positions at 2435 and 2114 cm~^ were consistent with 
D2O2 Jlannon et al.lll97lh . Assuming that the OD-stretching band was mainly derived from 



D2O2, the contributions of HDO to the band should be very small. 

As with the reaction OH + HD, two reaction channels were possible for the reactions 
of OD with HD: 

OD + HD ^ HDO + D, (20) 
OD + HD -> D2O + H. (21) 

For the reasons mentioned above, reaction (1201) was expected to be more efficient than 
reaction (ET]) . Figure El^a) shows an IR spectrum of the reaction product obtained after 
codeposition of D2O fragments with HD. Several new peaks were observed, most of which 
corresponded to HDO, as explained earlier. The column density of HDO was estimated from 
the OH-stretching band at 3400 cm~^ to be 3.3 x 10^^ molecules cm~^. 



- 15 - 



In contrast, the column density of D2O did not increase; rather, it decreased upon 
codeposition, as indicated by the decrease in peak area derived from the OD-bending band 
for D2O at 1200 cm~^ (Figure [H](b)). This result indicated that reaction (^T^ is not effective. 
Assuming that D2O was mainly produced through OD + OD D2O + O in the D2O- 
fragment deposition, the decrease in the D2O column density upon codeposition was mainly 
attributed to consumption of the corresponding OD by HD through reaction ( 120|) . This 
was based on the fact that reaction ( [T7|) (OD + D2 — > D2O + D) was not efficient, as 
demonstrated in Section 3.4. D2O2 was assumed to form via the reaction OD + OD, and 
it may be consumed in the presence of codeposited HD; however, variations in the D2O2 
column density could not be estimated because the OD bending band of D2O2 at 1000 cm~^ 
mixed with the bands corresponding to other species, such as O3 (Figure [8]^b)). 



3.6. Relative reaction efficiency 

We found that H2O and its isotopologues (HDO and D2O) form through surface reac- 
tions of OH or OD with H2, HD, or D2 at 10 K, and we determined the column densities of 
the reaction products in previous sections (Sections 3.1-3.5). The relative efficiency of each 
reaction was roughly estimated by comparing the column densities of the reaction products 
obtained under identical conditions (Table More accurate relative efficiencies, however, 
could not be obtained due to the reason shown below. As indicated in Section 2.1, the accu- 
rate fluxes of OH or OD were not measured experimentally, just assigned an upper limit of 
1.5 X 10^^ radicals cm^^ s^^, which is equal to the deposition rate of H2O or D2O when the 
microwave source is turned off. Although we verified by QMS and REMPI that almost all 
H2O and D2O molecules were fragmented when the plasma was turned on, unfortunately, the 
OH flux was not guaranteed to be identical to the OD flux, which makes the determination 
of the accurate relative efficiencies difficult in the present study. Nevertheless, the OH and 
OD fluxes did not vary signiflcantly from one another because the column densities of the 
reaction products, such as H2O or H2O2 in the H20-fragment deposition, were comparable to 
those of the isotopic counterparts, such as D2O and D2O2, in the D20-fragment deposition. 

The reactions listed in Table |2] are divided into three groups based on the efficiency 
of each reaction relative to reaction ( !T2|) : reactions with comparable efficiency (reactions 
( 1T5]) . (ITS]) , and ( 12U]) ). reactions one order of magnitude less efficient (reactions ( 1TB]1 and 
( 1T7|) ). and reactions of undetermined efficiency (reactions ( IT^ and ([21])). For convenience, 
hereafter, these three groups are denoted A, B, and C, respectively. Although the difference 
in the barrier height (AEa) among the reac tions listed in Table |2] is small [AE^ <~400 



K; iTalukdar et al.l Il996l : iNguyen et al.ll201l[ ). it is remarkable that a considerable isotope 
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effect was clearly observed in groups A and B. The reactions differed intrinsically depending 
on whether the reaction involved H- or D-atom abstraction, and the effective masses were 
comparable within a reaction group, as shown in Table |5J The effective masses for the 
reactions in group B (0.90) were almost twice as large as those in group A (0.47). The 
inverse correlation between efficiency and the effective mass of reactions is a typical feature 
of quantum tunneling because the tunneling rate decreases considerably as the effective mass 
of a reaction increases, as expected from Equation (fT3!l . 

The efficiencies of the reactions in group C can be determined from the increase in 
the column density of HDO or D2O in reaction f|T9|) or fl2T|) . respectively. The efficiency 
of reaction (|T9|) could potentially be estimated in a straightforward way by calculating the 
column density of HDO formed through the reaction. This was not possible, however, be- 
cause the OD-stretching band at 2435 cm~^, which was used to calculate the HDO column 
densities in the present study, included a significant contribution from D2O2, as explained 
above. The efficiency of reaction (12T]) could potentially be estimated from the increase in 
the D2O column density as a result of the reaction. The difference between the spectra of 
the codeposition and the D20-fragment deposition samples (Figure [S]^b)) shows, however, 
that the D2O column density did not increase; rather, it decreased with the codeposition of 
HD. We could not, therefore, calculate the increase in D2O column density through reaction 
f l2T|) . The efficiencies of the reactions in group C were thought to be comparable to those 
in group B because the effective masses are identical, i.e. the reactions in group C were less 
efficient than those in group A. A summary of the discussion in this section is schematically 
illustrated in Figure IHl which clearly shows that OH and OD preferentially abstract H atom 
from hydrogen molecules (H2 or HD), leading to H2O and HDO, respectively. 



ASTROPHYSICAL IMPLICATIONS 



The present study demonstrated that reactions between OH and H2 yiel d H^O on a solid 



surface at 10 K, despite a large activation barrier (2100 K) for the reaction (lAtkinson et al. 



20041 ). We assume that reaction ( [T2|) occurred in previ ous experiments on the formation of 



H2O through hydrogenation of O, O2, or O3 at 10 K (e.g. lMiyauchi et al.ll2008l : iMokrane et al. 



2OO9I ): however, this assumption was not extensively explored. This is the first experimental 
result to isolate reaction (IT^ and to sho w that reactions relat e d to H ^ actually yielded H2O. 
Although in a recent theoretical model, iDas fc Chakrabartil ( 1201 ll ) assumed that reaction 
f lT2|) does not occur on grain surfaces in dense clouds because of the activation barrier, the 
present experimental result strongly suggests that models should include reaction f lT2|) as 
the efficient route. 
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We propose that solid H2O in dense molecular clouds is mostly formed by the following 
three surface reactions: OH + H (reaction ([2])), OH + H2 (reaction f|T2|) ). and H2O2 + H 
(reaction ([S])). The main routes to the formation of H2O in dense clouds are schematically 
illustrated in Figure [101 OH would be formed through hydrogenation of O atoms and/or 
other reactions (e.g., reaction ([5])). It may be reasonable that H2O formation through OH 
+ H2 is more efficient than reactions OH + H in de nse clouds where H2 is more prevalent 
than H atoms among the various H-related species. iRomanzin et al.l (120111 ) also discussed 
the importance of reaction fll2p relative to reaction ([2]) for the formation of H2O in their 
experiments where energetic OH was produced by the hydrogenation of O3 at 25-50 K. As 
for H2O2, hydrogenation of O2 is probably a unique pathway toward the production of H2O2 
in dense clouds. O2 could be produced by the reaction of O + O on grain surfaces and/or 
supplied from the gas phase. The search for solid O2 in interstellar environments is ongoing 
and has not been successful to date. Its prese nce has been suggested in previous studies 



le.g. 



Boogert et al.ll2002l : [Goldsmith et al.ll201ll ). The recent detection of H2 O2 in the cloud 
core p Oph A suggested that reaction ([5]) may be plausible on grain surfaces (IBergman et al. 



2OIII ). The upper limit for the H2O2 ice abundance relative to H2O in several protostars 
and field stars was estimat ed based on labora tory IR spectra of H2O2/H2O mixed ice and 
was found to be 9% ± 4% (ISmith et al.l 1201 ll ). In addition, laboratory experiments clearly 
demonstrated that H2O form ation through O2 hydrogenation is effi cient under conditions of 



dense molecular clouds (e.g.. iMiyauchi et al.ll2008l : lOba et al.ll2009l ). 



In addition to the main routes mentioned above, various reactions were proposed that 
could lead to the formation of H2O and other species, such as H2O2 and OH (Table [T]). 
These reactions were assumed to make only small contributions to the observed abundance 
of solid H2O in dense clouds. For example, although hydrogenation of O3 was proposed 
as a n important route to t he formation of solid H^^O in dense clouds in previous theoreti 
cal (jTielens fc Hagenlll982l ) and experimental studies (IMokrane et al.ll2009l : iRomanzin et al. 
2OIII ). we suspect that the contribution of this reaction to H2O formation is small. O3 for- 
mation is not favored on grain surfaces because the mobility and flux of H atoms is higher 
than the corresponding properties of O atoms, which prohibits formation of O3 by consuming 
O atoms (reaction ([T])) and/or O2 (reaction ([3])) prior to O3 format i on. O ur conclusion is 
consistent with the theoretical model proposed by lCuppen fc HerbstI (120071 ). who estimated 
the efficiencies of reactions ([2]), ([5]), and (fT2!) in dense clouds at 10 K to be 6%, 17%, and 
77%, respectively. 

The formation of solid HDO was experimentally demonstrated to occur in the present 
study. Although HDO in the sohd state has been sought after, its detection has not been 
successful to date. An upper limit for t he solid HDO/H2O ratio in some sources is available 
(jPartois et al.ll2003l : iParise et al.ll2003l ). In the gas phase, the HDO/H2O ratios in the low- 
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mass protostars IRAS 16239-242 2 and NGC1333-IRAS2A were found to be 0.03 flParise et al. 
2005h and 0.07 Jliu et al.lboill ). respectively. 



Kristensen et al.l ( 120111 ) discussed the formation of H2O, HDO, and D2O through reac- 



tions of OH with H2, HD, and D2, respectively, on ice-covered dust grains in dense molecular 
clouds. They assumed that the reaction OH + HD is barrierless, and proposed that the 
H2O and HDO form statistically in a ratio of 2:1, respectively, from the intermediate H2DO. 
Based on these assumptions, this group claimed that the astronomically observed HDO/H2O 
ratio could be explained according to their model; however, the reactions of OH/OD with 
H2/HD/D2 do have a bar rier, as demonstrated in previous experime ntal and theoretical stud- 



ies in the gas phase (e.g.. iTalukdar et al.lll996l : iNguyen et al.ll201ll ). as well as in the present 



study. Moreover, the present study demonstrated that the formation of HDO through the 
reactions of OH with HD were much less favored than the formation of H2O, which is incon- 
sistent with their assumptions. Our results are i nconsistent with the theor etical predictions 
for the formation of H2O and its isotopologues by lCazaux et al.l (120101 . 12OIII ). who considered 
branching ratios of 1 to 1 between reactions (fTSi) and (fT9!) or between reactions (!20l) and 

The present experimental results demonstrate that a crucial factor for constraining 
the ratio of the HDO and H2O abundances formed by reactions listed in Table [2] is the 
relative abundance of OH and OD in dense clouds. The following three reasons support 
this conclusion: first, the efficiency of H2O formation through the reactions of OH with H2 
(reaction (IT^ ) and HD (reaction (ITSl) ) is comparable to the efficiency of HDO formation 
through reactions of OD with H2 (reaction (fT5l) ) or HD (reaction (120|) ). Second, the four 
reactions (i.e., H-atom abstraction reactions) are much more efficient than the other four 
reactions (i.e., D-atom abstraction reactions) (Table [2], Figure [9]). Third, on ce solid H-^O is 



formed, its deuterati on does not occur t hrough H-D exchange with D atoms (INagaoka et al. 



20051 ) or with D2O (IGalvez et al.l 1201 ll ) at 10 K, even if D2O ice forms successfully. Note 
that this prediction is only applicable to HDO and H2O formation through reactions listed 
in Table H 

If OH and OD are mainly produced by the reactions of O atoms with H and D atoms, 
respectively, the relative abundance of H and D atoms would be important to constrain the 
relative abundance of OH and OD. Accord ing to the deuterium fractionation models in the 
gas phase proposed by [Roberts et al.l (120021 ). the atomic D/H ratio is low (< 10~^) during the 
early stages (< 10^ years) of molecular cloud evolution. After 10^ years, the ratio increases 
dramatically, reaching a value of ~0.06 after 10^ years, which is c omparable to the observed 
HDO/H2O ratio shown above (jParise et al.ll2005l : iLiu et al.ll201ll ). If large regions of an ice 
mantle are produced during the later stages of molecular cloud evolution, the HDO/H2O 
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ratio is expected to reach its highest value for the D/H atomic ratio, ~0.06 ( iRoberts et al. 



20021 ). If not, other routes may contribute to the observed HDO/H2O ratio. 



The above scenario is based on a n assumption that O H (OD) is produced through the 



reaction O + H (D) on grain surfaces. [Roberts et al.l (120021 ) estimated the OD/OH ratios in 



their gas-phase model to be 0.093-0.358, higher than the atomic D/H ratio in the same model. 
The participation of gas-phase OD and OH in deuterium fractionation in water formation 
via surface reactions leads to D enrichment in ice, although the degree of enrichment will 
depend on the flux of the gas-phase OD and OH onto grains. 

Further studies are required to inclusively discuss the formation of solid H2O and HDO 
in dense molecular clouds. In future studies, care should be taken when considering chem- 
ical reactions, particularly quantum tunneling reactions, with comparable activation barri- 
ers because the reaction efficiency can differ dramatically among such reactions, as clearly 
demonstrated in the present study. 
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Fig. 1. — (a) (2 + 1) REMPI spectrum of OH emitted from the radical source through 
the Al pipe at room temperature and (b) the best-fit simulated spectrum for Trot = 300 
K. The capital numbers 1-6 indicate rotational levels of Q, R, and S branches associated 
with the X^U ground state of OH. The subscript numbers 1 and 2 represent the lower or 
upper members of the rotational level in the X^U ground state, respec tively, because spin- 



orbit coupling causes the rotational levels to split into two components (Ide Beer et al.lll991 



Dieke fc Crosswhitelll962[ ). 
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Fig. 2. — FTIR spectra after (a) deposition of H2O fragments and (b) codeposition of H2O 
fragments with H2. (c) The spectral differences relative to the H2O- fragment deposition 
sample. The peak with the asterisk (*) was used for quantification purposes. The peaks 
at ~2300 cm~^ and ~2700 cm^^ were derived from the background CO2 and the inherent 
noises caused by vibration of the He refrigerator, respectively. 
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Fig. 3. — FTIR spectra after (a) deposition of D2O fragments and (b) codeposition of D2O 
fragments with H2. The peak with the asterisk (*) was used for quantification purposes. 
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Fig. 4. — Integrated TPD yield of HDO {m/z — 19) as a function of the IR peak area of 
the 1450 cm^^ band, which appeared upon the codeposition of D2O fragments with H2. The 
data points show the averaged values with statistical errors for four different durations. A 
dashed line represents the best-fitting straight line through the points. 
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Fig. 5. — FTIR spectrum after codeposition of H2O fragments with D2. The peak with the 
asterisk (*) was used for quantification purposes. The inset shows the OD-stretching band 
of HDO at 2475 cm-^ 
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Fig. 6. — FTIR spectrum after codeposition of D2O fragments and D2. The peak with the 
asterisk (*) was used for quantification purposes. 
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Fig. 7. — (a) FTIR spectrum after codeposition of H2O fragments witli HD and (b) tlie 
spectral differences relative to the H20-fragment deposition sample. The peak with the 
asterisk (*) was used for quantification purposes. The peaks at 2300 cm~^ were derived from 
the background CO2. 
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(a) D^O fragments + HD 
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Fig. 8. — (a) FTIR spectrum after tlie codeposition of D2O fragments witli HD and (b) 
the spectral differences relative to the D20-fragment deposition sample. The peak with the 
asterisk (*) was used for quantification purposes. 
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Fig. 9. — Schematic illustration showing the reaction network for the formation of solid H2O 
and its isotopologues (HDO and D2O) via the reactions listed in TableO The broad or narrow 
lines represent reactions with relative efficiencies of about 1 or 0.1, respectively, compared 
to reaction f|T2l) (Table [2]). The dashed lines indicate reactions with efficiencies that could 
not be estimated experimentally. The numbers correspond to the reaction numbers given in 
the text and Table [2l This figure clearly shows that OH and OD preferentially abstract H 
atom from hydrogen molecules (H2 or HD), leading to H2O and HDO, respectively. 
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Fig. 10. — Main routes to the formation of solid H2O in dense molecular clouds: liydro- 
genation of OH (reaction (jj])), the reaction of OH with H2 (reaction (JT2l) ). and the reaction 
of H2O2 with H atom (reaction ([5])), which are denoted by solid lines. Other reactions are 
denoted by dotted lines. The title reaction OH + H2 is denoted by a thick line. The numbers 
in parentheses correspond to the reaction numbers given in the text and Table [TJ 
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Table 1. Surface reactions proposed for the formation of solid H2O and the related species 



Number 


Reaction 


Ea (K) 


AH°{kJ mol-i) 


Reference 


1 


+ H ^ OH 





-429 




2 


OH + H ^ H2O 





499 




3 


O2 + H ^ HO2 


~0 


-203 


a 


4 


HO2 + H ^ H2O2 





-354 


b 


5 


H2O2 + H -)■ H2O + OH 


2000 


-285 


c 


6 


OH + OH H2O2 





-211 


d 


7 


OH + OH -> H2O + 





-71 


d 


8 


HO2 + H ^- H2OO* H2O + 0(iD) 





-lO'' 


b 


9 


HO2 + H ^ H2O..O ^ H2O + 0(3p) 


9300 


-242 


b 


10 


HO2 + H ^ H2 + 02(3E -) 


750 


-229 


b 


11 


O3 + H ^ OH + O2 


500 


-321 


e 


12 


OH + H2 -)■ H2O + H 


2100 


-62 


d 


22 


+ H2 OH + H 


3160 


+8 


f 


23 


HO2 + H2 H2O2 + H 


13100 


+58 


g 



^Walch et al. (1988); Sellevag et al. (2008). 

'•Mousavipour & Saheb (2007). 

'^Koussa et al. (2006). 

d Atkinson et al. (2004). 

^Kcyscr (1979). 

fBaulch et al. (1992). 

STsang & Hampson (1986). 

''This rea<;tion as a whole is exothermic; however, the latter part of the reaction (H2OO* — ^ 
H2O + 0(iD)) is endothermic by 151 kj mol"!. 

Note. — AH° denotes the experimental enthalpy change for the gas-phase reactions in kJ 
mol^^, where 1 kJ mol~^ corresponds to 120 K. Values of Ea and AH° are derived from gas- 
phase reactions. 
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Table 2. Summary of the Activation Barrier (K), the Effective Mass, and the Measured 
Efficiency for Reactions Off + ff2 and the fsotopologues. 



Number 




Reaction 




Section^ 


b 


Ea (K) 

c 


d 


Atom 
Abstracted 


Effective 
Mass" 


Relative 
Efficiency' 


12 


OH - 


h Ha HaO 4- 


H 


3.1 


2000 


2100 


2935 


H 


0.47 


1 


15 


OD - 


h H2 ^ HDO H 


h H 


3.2 






2789 


H 


0.48 


~1 


16 


OH - 


h D2 -s> HDO H 


h D 


3.3 


2456 




3026 


D 


0.90 


~0.1 


17 


OD - 


h D2 D2O 4 


D 


3.4 






2870 


D 


0.91 


~0.1 


18 


OH - 


h HD -5- H2O H 


h D 


3.5 


2130 




2855 


H 


0.48 


~1 


19 


OH - 


h HD -s> HDO 


+ H 


3.5 


2130 




3051 


D 


0.90 


n.d.s 


20 


OD - 


h HD HDO 


+ D 


3.5 






2703 


H 


0.48 


-1 


21 


OD - 


h HD D2O - 


h H 


3.5 






2900 


D 


0.90 


n.d.s 



^A corresponding section in the present paper. 

Experimentally determined values jTalukdar et al.lll996l 'l. 
'^Compilation of literature values l l Atkinson et al.ll2004h . 

Computed reaction barrier heights including harmonic zero-point vibrational energy llNguven et al.ll201 j) . 
"Calculated by using Equation 1141 1 l |johnstonlll966l ). 

'Estimated by comparison of the yield of the product (H2O, HDO, or D2O) for each reaction with that of H2O for 
reaction II12I I. 

^Not determined due to reasons described in the text. 



